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Abstract — We develop distributed algorithms to allocate re- 
sources in multi-hop wireless networks with the aim of mini- 
mizing total cost. In order to observe the fundamental duplexing 
constraint that co-located transmitters and receivers cannot oper- 
ate simultaneously on the same frequency band, we first devise a 
spectrum allocation scheme that divides the whole spectrum into 
multiple sub-bands and activates conflict-free links on each sub- 
band. We show that the minimum number of required sub-bands 
grows asymptotically at a logarithmic rate with the chromatic 
number of network connectivity graph. A simple distributed and 
asynchronous algorithm is developed to feasibly activate links 
on the available sub-bands. Given a feasible spectrum allocation, 
we then design node-based distributed algorithms for optimally 
controlling the transmission powers on active links for each sub- 
band, jointly with traffic routes and user input rates in response 
to channel states and traffic demands. We show that under 
specified conditions, the algorithms asymptotically converge to 
the optimal operating point. 



I. Introduction 

While offering the potential for ubiquitous and untethered 
communications, wireless networks typically demand more 
sophisticated resource management than wireline networks. 
Optimal resource allocation in large-scale wireless networks 
involves joint optimization across multiple layers as well 
as distributed implementation across network nodes. In this 
paper, we develop distributed algorithms which jointly allocate 
frequency spectrum, ttansmission powers, traffic input rates, 
and traffic routes on a node-by-node basis in order to mini- 
mize total cost in an interference-limited multi-hop wireless 
network. 

While joint optimization involving power control, conges- 
tion control, and routing has been studied previously [1]— [3], 
a common shortcoming of the previous work is the failure 
to account for the constraint that a node cannot transmit 
and receive simultaneously on the same frequency band. For 
while some types of multi-user interference can be ameliorated 
by using advanced coding techniques, interference between 
a transmitter and a co-located receiver is very difficult to 
suppress due to the transmitted power being many orders 
of magnitude higher than the received power at the same 

'This research is supported in part by NSF grant CNS-0626882 and AFOSR 
grant FA9550-06-1-0135. 



node0 The resulting constraint that a transmitter cannot be 
simultaneously active with a co-located receiver [4] on the 
same frequency band is referred to as the duplexing constraint. 
In practical wireless networks, the duplexing constraint ap- 
pears to be quite fundamental, and thus must be observed by 
resource management schemes. 

In this paper, we develop distributed resource allocation 
algorithms for wireless networks in accordance with the 
duplexing constraint. To accomplish this, we first devise a 
new spectrum allocation scheme which divides the spectrum 
into a sufficient number of frequency bands and activates 
co-located transmitters and receivers on different bands. We 
show that the minimum number of sub-bands needed for 
resolving duplexing conflicts is asymptotically logarithmic in 
the chromatic numbei0 of the network connectivity graph. We 
provide a simple algorithm that feasibly assigns frequency 
bands to links in a distributed and asynchronous manner with 
low control overhead. 

Given a conflict-free spectrum allocation, we then design 
a set of node-based distributed gradient projection algorithms 
that iteratively adjust transmission powers, traffic input rates, 
and traffic routes according to channel conditions and traffic 
demands, in order to minimize total network cost. The power 
control and routing algorithms we develop are frequency 
selective, in that for each link, the power conttol algorithm 
adjusts the transmission power on each of the link's active sub- 
bands. The routing algorithm involves both inter-node routing, 
which specifies the allocation of incoming traffic at each node 
to its outgoing links, and intra-node routing, which specifies 
the allocation of the total traffic on a given link across its 
active sub-bands. Finally, we show that congestion control 
can be naturally incorporated by considering an equivalent 
routing problem on a virtual overflow link. We show that under 
specified conditions, the iterative algorithms converge to the 
optimal operating point from any initial condition. 

theoretically a node is able to subtract the transmission signals generated 
by itself from the received signals so that the self-interference can be perfectly 
cancelled out. In real decoders, however, the received signals are practically 
irrecoverable in the face of the overwhelming transmitted signals. 

3 The chromatic number of a graph is the minimum number of colors with 
which the vertices of the graph can be colored such that adjacent vertices 
have different colors. 
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The rest of the paper is organized as follows. In Sec- 
tion in we introduce the network model, discuss the duplexing 
constraint, and formulate the spectrum allocation problem. 
In Section [HI] we find the minimum number of sub-bands 
required by a feasible spectrum allocation. A distributed and 
asynchronous sub-band allocation algorithm is developed in 
Section [TV] In Section [V] we formulate the cross-layer op- 
timization problem for networks operating on multiple sub- 
bands, and derive the conditions satisfied by the optimal 
configuration. In Section [VI] we present node-based gradi- 
ent projection algorithms to jointly optimize power control, 
routing, and congestion control based on the outcome of the 
spectrum allocation, as well as the channel conditions and 
traffic demands of the network. 

II. Network Model and Spectrum Allocation 

Let the wireless network be modelled by a directed and 
connected graph Q = (Af, C), where Q is referred to as the 
connectivity graph of the network. A node i e M represents a 
wireless transceiver, and G C represents a unidirectional 
wireless link from node i to j. We assume Q is link-symmetric, 
i.e., if G C, then (j, i) G C, and vice versa. Two nodes i 
and j are neighbors if G C (or equivalently G C). 
Let Oi = {j : G £} be the set of i's neighbors. 

A. Duplexing Constrained Interference Model 

We focus on a network model that prohibits any node 
from simultaneously transmitting and receiving on the same 
frequency band. That is, we impose the duplexing constraint 
on every band. Such a constraint is less stringent but more 
fundamental than the extensively studied primary interference 
constraint [5]-[7], where any node can transmit or receive 
(at any given time and on any given band) on at most 
one active link, and the secondary interference constraint, 
which further prohibits any node from transmitting when there 
is a neighbor receiving from another node [8]. A general 
approach to studying a broad class of interference constraints 
was presented in [9]. Indeed, duplexing constraints cannot be 
bypassed (at least currently) by using sophisticated coding 
methods, and must be observed by practical network man- 
agement techniques. 

Traditional network management techniques which aim to 
resolve various types of interference among wireless links have 
concentrated on scheduling in time [10]— [17], where at any 
given time, only mutually-non-interfering links are activated. 
Scheduling, however, usually requires centralized controllers 
and involves high communication and computational complex- 
ity [7], [18]. Simplified distributed scheduling policies have 
been proposed [5], [6], [8], [19], [20] for various purposes. 
In general, however, the reduced implementation complexity 
comes at the expense of performance [21]. 

The difficulty in finding interference-free schedules in time 
leads us to seek an alternative solution. A natural approach 
is to consider network management in the domain of fre- 
quency instead of or in addition to time [22]-[25]. Because 
communication on different frequency bands are practically 
non-interfering, one can think of simultaneously applying 



different link activation sets on non-overlapping frequency 
bands within the assigned spectrum. In this scheme, nodes 
transmit on certain bands while receiving on other bands to 
avoid duplexing interference. The spectrum allocation tech- 
nique has an important advantage over scheduling in time: 
once a feasible spectrum allocation is established, the network 
can function relatively statically in that mode, not having to 
switch to another mode unless network itself changes substan- 
tially. While spectrum allocation problems have been proposed 
and studied in the interference graph induced by particular 
interference constraints, the solution using existing vertex- 
coloring methods [22], [23] is cumbersome and its complexity 
scales poorly with the size of the network. Moreover, the 
optimization of frequency assignment techniques has not been 
thoroughly investigated. In particular, the number of available 
frequency bands is often arbitrarily set and frequency bands 
are assigned to links in a heuristic manner [24], [25]. 

In this work, we adopt the spectrum allocation approach 
to resolve the fundamental duplexing conflicts for general 
wireless networks. In particular, we investigate two central 
questions: (1) what is the minimum number of frequency 
bands with which all co-located transmitters and receivers can 
be simultaneously activated subject to duplexing constraints, 
and (2) how can one efficiently find a feasible frequency 
assignment when there are enough frequency bands? We 
provide an exact analytical answer to question (1) and develop 
a distributed asynchronous algorithm which solves problem 
(2). Our analysis is based only the network connectivity graph. 
This approach requires much less storage and computation 
overhead than alternative methods that utilize the interference 
graph. 

B. Spectrum Division and Sub-band Allocation 

The duplexing constraint permits only a subset of the links 
to be activated simultaneously on each frequency band. To 
activate all links at the same time, it is necessary to divide the 
spectrum into several sub-bands and activate different subsets 
of conflict-free links on different sub-bands. 

Suppose the network occupies a contiguous spectrum which 
can be partitioned into a number, say Q, of sub-bands, each 
of which covers a contiguous segment of the whole spectrum. 
Let the collection of the sub-bands be denoted by Q. With a 
specific spectrum division in place, each link can be active on 
one or more of the sub-bands. If link uses sub-band q, 
we say is an active link on q, and q is an active sub- 

band of (i, j). Denote the subset of links that are active on a 
sub-band q by C q , and the set of active sub-bands of by 
Qij. A spectrum allocation is given by the collection {£ q } q eQ 
(or equivalently {Qij}(i,j)ec)- Note that finding a spectrum 
allocation involves two steps: spectrum division, which decides 
how many sub-bands the whole spectrum is divided into, and 
sub-band allocation, which determines which links are active 
on which sub-bands. We will address these two issues in 
Sections [III] and IIVI respectively. 

A spectrum allocation is feasible if (i) for all G C, 

Qij 0, and (ii) for all q G Q, C q satisfies the duplexing 
constraint. Thus, any node's outgoing and incoming links 
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cannot be both active on the same sub-band. However, it is 
feasible for a node to have multiple active outgoing links or 
multiple active incoming links on a single sub-band. 

C. Interference Graph and Number of Sub-bands 

In previous studies on frequency assignment techniques, 
the minimum number of frequency bands is found from 
the interference graph induced by the specific interference 
constraints [22], [23]. For the duplexing constraints, the in- 
terference graph Q is constructed as follows. Let the vertices 
of the interference graph Q = (V, £) be the links of the 
network connectivity graph Q = (Af, C), i.e., V = C. In Q, 
an edge exists between two vertices (links in Q) if one link's 
transmitter is the other link's receiver. The interference graph, 
unlike the connectivity graph, is undirected^ It is easy to see 
that a feasible spectrum allocation on Q exists if and only if 
the number of available sub-bands is greater than or equal to 
the chromatic number x(G) of the interference graph Q. We 
illustrate the interference graph approach in Figure Q] using a 
complete four-node network. On the left is the connectivity 
graph Q, whose induced interference graph Q is shown in 
the middle. A minimal vertex-coloring on Q using four colors 
{R, G, B, Y} is depicted. The right graph represents the link- 
coloring on Q inferred by the minimal vertex-coloring on Q. 
The color(s) associated with each node are those assigned to 
the outgoing links of that node. The reason for using such a 
representation will be explained shortly. 

Although the interference graph approach can provide an 
answer to our problem, it has major shortcomings. Note that 
the size of the interference graph |V| = \C\ is on the order 
of | A/] 2 . To compute the chromatic number of Q, Q has to 
be constructed and stored at a central controller which then 
computes x(G) by finding a minimal vertex-coloring, which 
is itself an NP-complete problem [26]. For these reasons, the 
interference graph approach is not tractable for medium-and- 
large-scale networks. This motivates us to seek an alternative 
method that can find the minimum number of sub-bands 
directly from the connectivity graph Q. 

III. Spectrum Division 

In this section, we investigate the minimum number of sub- 
bands which yields a feasible spectrum allocation. Our analysis 
will be based entirely on the network connectivity graph Q. 
For convenience of exposition, we first transform the problem 
into an equivalent graph-theoretic link-coloring problem. 

A. Link-Coloring Problem 

For the moment, we leave the total number of sub-bands 
undetermined. Let each sub-band be identified by a unique 
color. We represent a spectrum allocation by a color assign- 
ment to links. A feasible spectrum allocation corresponds to 
a color assignment such that (i) all links are assigned with at 
least one color and (ii) for any node, no incoming link has a 

4 In order to avoid confusion, we refer to the connectivity graph as having 
nodes J\f and links C, and the interference graph as having vertices V and 
edges £. 



common color with any one of its outgoing links. We will find 
the minimum number of colors required for such a coloring. 
For any color assignment, denote the set of colors used by 
node i's outgoing links by OCi, and the set of colors used by 
its incoming links by IC L . Since the network is connected and 
link-symmetric, every node has at least one outgoing and one 
incoming link, implying that OCi and ICi are both nonempty 
for all i G Af. A color assignment is feasible if and only 
if OCi H ICi = for all i G AfB Note that the feasibility 
of a coloring scheme can be checked using the node color 
sets {OCi : ICi]i(zj^, regardless of the detailed assignment 
of colors to links. In fact, we will show that studying the 
node outgoing color sets {OCi} alone suffices for the link 
coloring problem. In the following, we say a color assignment 
is consistent with {OCi} if {OCi} results from the color 
assignment. 

Lemma 1: Given a graph Q = {Af, C} and the nodes' 
outgoing color sets {OCi}i^, there exists a consistent and 
feasible color assignment if and only if 

(i) Od\OC 3 ^ 0, for all G £; 

(ii) U iGOl OC\OCj = OCi, for all i e Af. 

Proof: We first prove the necessity part. Suppose first that (i) 
is violated, i.e., there exists a link such that OCi Q OCj, 
then must have a color in OCj, meaning that some 

outgoing link of j has the same color as which is 

infeasible. Next suppose that (ii) does not hold for some node 
i, i.e., there exists a color q G OCi such that q G OCj of 
every j G OCi. Since q is used on some outgoing link of i, 
say (i, k), (i, k) will have a color belonging to OCk, violating 
the duplexing constraint at node k. 

The conditions are sufficient because given {OCi}i^ 
which satisfy (i)-(ii), a feasible coloring can be constructed 
by assigning each link the color(s) in OCf\OCj. By 

(ii), this coloring is also consistent with {OCi}i^ . □ 

In the following, we use the shorthand notation OCi ^ OC j 
to represent OC l \OC j ^ and OC 3 \Od ^ 0. A collection 
of node outgoing color sets {OCi}i^M is said to be feasible 
if it satisfies the conditions in Lemma [T] 

We now state the problem of finding the minimum number 
of sub-bands in terms of the minimal coloring problem as 
follows. Given a graph Q — (Af,C), 



minimize 



subject to 



U Od 
{OCi}i e j^f feasible. 



(1) 



Denote the minimum number of colors by Qg. To solve for 
Qg, we first relax the constraints of problem (Q~|) by dropping 
the second condition in Lemma Q] and consider 



minimize 



subject to 



BtocBj, V(i,j)6£, 



(2) 



Note that we allow one link to transmit on multiple sub-bands. Hence, we 
do not require \Od\ < \Oi\ or \ld\ <\Gi\. 
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where Bi, i G Af are any non-empty sets. The optimal solution 
Qg to (O must be less than or equal to Qg. We will find a 
lower bound on Qg and a matching upper bound on Qg, thus 
uniquely determining Qg. For lower bounding Qg, the next 
observation is useful. The proof is deferred to Appendix lAl 

Lemma 2: Let B\ , • ■ • , _Bfc be distinct subsets of a Q-set 
all of which have the same cardinality g. If g > |_Q/2J, then 
there exist distinct subsets G\ , ■ ■ ■ ,Ck such that Cj C Bj and 
|Cj| = g - 1, i = 1, • • • , k. If g < \Q/2\, then there exist 
distinct subsets C\,--- ,Ck such that Bj C Cj and |Cj| = 
.9+ 1, i = !=••• 

Using Lemma [2] and the fact that any two distinct subsets 
Bi, Bj with the same cardinality satisfy Bi oo Bj, we can 
show the following (see the proof in Appendix IB!. 

Lemma 3: There exists an optimal solution {Bf}i e ^ for (O 
such that all Bi, i 6 Af, have the same cardinality. 

Using Lemma [3] we can find a lower bound on Qg as 
follows. 

Lemma 4: Given a graph Q = (Af, C) with chromatic 
number x(G)' the solution Q*g of the problem (O is greater 
than or equal to Q(x(G)), where the function Q : Z + i— ► Z + 
is defined as 



Q(N) = min <^ q E Z+ 



q 

L9/2J 



> TV 



(3) 



Proof: By Lemma [3] we can without loss of optimality 
consider a solution {Bf}i e ^ of (O such that |£^| = \Bj\, 
for all i,j G jV. Therefore, ^ Bj whenever Bi ^ 
Suppose there are K distinct elements in {Bi}i £ j^. Because 
lUi^i! = Q*g> by Sperner's theorem [27], if < ( LQ ^ 2J ). 
Moreover, since {Bi} satisfies the constraint in (O, we must 
have K > x(^), where x(^) i s m e chromatic number of Q. 
Therefore, Q* g > Q{ X {Q))- □ 

Lemma |4] indicates that at least Q(x(Q)) colors are needed 
to construct a feasible collection of node color sets even with 
the condition (ii) in Lemma[T]removed. However, by assigning 
node color sets according to a minimal node-labelling scheme 
of Q, Q(x(G)) colors are sufficient to yield a feasible configu- 



ration satisfying both conditions in LemmaQ] We can therefore 
conclude the following. 

Lemma 5: Given a graph Q = (Af, C) with chromatic 
number x(S), the solution Qg to the problem (Q3 is less than 
or equal to Q(x(Q))- 

Proof: For graph Q, Q(x(Q)) colors are sufficient because 
they can generate at least x(^) color subsets {Bi} with 
equal cardinality \Q(x(Q)) /2J, any two of which satisfy 
Bi no Bj. Associate any x(Q) of the color subsets with x(Q) 
labels. Assign the x{Q) color subsets to nodes according to a 
minimal node-labelling scheme, which requires exactly x{Q) 
labels. We then have a feasible collection of node outgoing 
color sets that satisfies the two conditions in Lemma [Trl □ 

We have from Lemmas [4] and that 

Q(x(Q)) <Q* g <Qg <Q(x(Q))- 

Therefore, all inequalities must hold with equality. 

Proposition 1: Let the chromatic number of graph Q be 
x(Q)- Then, the solution of (HJ is Qg = Q(x(Q))^ where 
the function Q(-) is given by OJ. 

It is easy to see that Q(-) is nondecreasing. The values of 
Q(N) for N = 1, • ■ • , 20 are plotted in Figure |2] For large 
N, it follows from Stirling's approximation [28] 



q 

VL?/2J 



exp 



that 



Q(N) = e(logN). 

Therefore by Proposition Q] Qg grows at a logarithmic rate 
with x(Q)- 

IV. Distributed Sub-band Allocation 
A. Spectrum Division Using Approximation of x(Q) 

We have found that the minimum number of sub-bands 
depends on the chromatic number x(S)- The problem of 

6 An arbitrary color subset assignment according to a node-labelling solution 
generally results in {OC;}igj\f satisfying only (i) in Lemma [T] Should (ii) 
not hold for some i, reset OCi := Ujgo OCi\OCj . With this modification, 
both (i) and (ii) in Lemma [T] hold and no extra color is needed. 
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Fig. 2. The values of Q(N), N = !,-■■ ,20. 



finding x(G)< however, is NP-complete [26] in general. On 
the other hand, any upper bound on x{Q) gives us a sufficient 
number of sub-bands. One such well-known upper bound is 
x(Q) < A(Q) + 1, where A(Q) is the maximum degree of 
any node in Q. In the following, we assume that there are at 
least Q(A(Q) + 1) sub-bands available. 

The maximum degree A(Q) is straightforward to determine 
(in a distributed manner), and we assume that it is known to 
all nodes a priori. Notice that Q(A(Q) + 1) is not too far 
from Q(x(S)) in typical networks^ since the Q(-) function 
is piecewise flat and grows approximately as a log function. 
It is worth comparing the upper bound Q(A(Q) + 1) with 
A(Q) + 1, an upper bound of x(G), the chromatic number of 
the interference graph. Recall that vertices of Q are (directed) 
links of Q. Any vertex («, j) of Q has neighbors {(fc, I) € C : 
k = j or I = i, but not both}. Hence, by the assumption that 
Q is link-symmetric, the degree of in Q is given by 

A (lj) = \Ii\ + \Oj\ — l = Ai + Aj - 1. 

Here, A, and Aj are the degrees of nodes i and j in Q. 
Consequently, the maximum degree of Q is 

A(g) = max Aj + A, - 1. 

(<,j)ec 

On the other hand, A{Q) = maxi^A^. Therefore, it is 
always true that 

A(G) > A(G). 

The equality holds only when the network topology is ex- 
tremely asymmetric, e.g. the star networks. Typically, A(Q) 
can be almost twice as large as A(Q). Moreover, the Q(-) 
function is piecewise flat and scales as log 2 ( ). The difference 
of the two upper bounds A(Q) + 1 — Q(A(Q) + 1) typically is 
substantial. Therefore, the connectivity graph approach devel- 
oped in this work not only is more simple and straightforward 
but also leads to a tighter upper bound than the conventional 
interference graph approach. 

B. Distributed Sub-band Allocation Algorithm 

With the spectrum division accomplished, we now devise 
an algorithm with which feasible {OCi} can be determined 

7 There are graphs, e.g. complete graphs, for which x(G) < A(CJ) + 1 
holds with equality. 



in a distributed manner and a feasible spectrum allocation can 
be found. The algorithm applies to all connected and link- 
symmetric graphs. Assume that Q sub-bands are available. 
For expositional purposes, we keep an unprocessed node set 
U, which initially contains all the nodes. The Distributed Sub- 
band allocation (DSA) algorithm is iterated as follows. 

Step 1. Initially U = N. Arbitrarily select a node i from U 
and set U := U\{i}. Arbitrarily choose [Q/2J sub-bands to 
form OCi, mark i as "processed", and go to Step 2. 

Step 2. Select an arbitrary node i from U such that i has 
at least one processed neighbor (under the assumption that 
the network is connected, there is always such a node in hi 
after the first node is processed) and set U := U\{i}. Node i 
finds an OCi with [Q /2J sub-bands which is different from all 
OCj of its processed neighbors j (it can be shown that such 
an OCi always exists). Moreover, to maximally avoid potential 
channel interference, i selects sub-bands in increasing order of 
their number of occurrence (gin {OCj} j& 0i \u (ties are broken 
arbitrarily)]^ Mark i as "processed". If U ^ 0, repeat Step 2; 
otherwise go to Step 3. 

Step 3. Each node i allocates the sub-band(s) in Qjj := 
OCi\OCj to each outgoing link (i, j). The algorithm termi- 
nates. 

Proposition 2: The set {Qij}^ ,j)eC generated by the DSA 
algorithm induces a feasible spectrum allocation {£ 9 } g6 g. 

Proof: First we show that in Step 2 of the algorithm, an 
appropriate OCi can always be found. Because the whole set 
has Q > Q(A(Q) + 1) elements, there are at least A(Q) + 1 
distinct subsets with cardinality [Q/2\. However, i has at most 
A(Q) (processed) neighbors. Hence, there always exists at 
least one candidate for OCi. 

Now we show that the {OCi} obtained up to Step 3 are 
feasible. They satisfy (i) in Lemma Q] due to the rule of 
successively selecting OCi in Step 2. They also satisfy (ii) in 
Lemma [T| as proved next. Suppose on the contrary that there 
exists i e M and q e OCi such that q <E OCj for all j <E O,;. 
If i is not the first processed node, the fact that q G OCi 
implies that by the time OCi is being determined, each of the 
Q ~ LQ/2J sub-bands not included in OCi must appear in 
every previously processed neighbor j's OCj. Taking also q 
into account, we can deduce that all those OCj have at least 
LQ/2J + 1 elements, which is a contradiction. Even if i is 
the first processed node, the next processed node j must be 
i's neighbor and OCj U OCi = due to the algorithm's rule. 
Hence, the hypothesis cannot be true in either case. 

With the feasibility of the {OCi} established, it is easy to 
verify that the {Qij}(i.j)£c generated in Step 3 induces a 
feasible spectrum allocation {Cq} q eQ- D 

In practice, the DSA algorithm can be implemented in a 
distributed fashion by nodes in the network. Specifically, after 
any node i has arbitrarily set its OCi at the beginning, any 

8 The number of occurrences can be zero. 

9 For instance, if the whole set of sub-bands is {R, G, B, Y}, and node i 
finds its two neighbors have chosen {R, G} and {R, Y}, it will accordingly 
choose either {B, G} or {B, Y}. 
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other node j can determine its own OCj as long as it has at 
least one processed neighbor and no other neighbor is being 
processed at the same instant. In other words, it is possible to 
have two non-adjacent nodes configuring their outgoing sub- 
bands at the same time. Since the action of a node depends 
only on its neighbors, node operations need not be globally 
coordinated or synchronized across the network. All that is 
required is an initialization phase that designates the first node 
to be processed. Thus, the sub-band allocation algorithm is 
distributed and asynchronous in nature. 

Finally, we note that the DSA algorithm is robust to dy- 
namic node failure and addition. Suppose that after a feasible 
spectrum allocation is established by the DSA algorithm, a 
node leaves the network. In this case, the configuration for 
the remaining nodes and links is still feasible. If a new node i 
is to join the network, it can run Step 2 of the DSA algorithm 
to determine its own OCi based on its neighbors' OCj,j S Oi. 
Then, active sub-bands for each link incident to i can be 
allocated as in Step 3. As long as i connects to at most A(Q) 
nodes in the old network Q, the algorithm yields a feasible 
spectrum allocation for node i, with all other nodes' allocation 
unaffected. 

V. Cross-layer Optimization in Multi-radio 
Multi-hop Wireless Networks 

The algorithm in Section[lV]generates a spectrum allocation 
where the active links on every sub-band satisfy the duplexing 
constraints. Given this feasible spectrum allocation, we now 
develop corresponding joint power control, routing, and con- 
gestion control algorithms which minimize total network cost, 
which will be specified later, given channel conditions and 
traffic demands. 

A. Interference Limited Transmissions and Node Power Con- 
straints 

We consider networks where messages on each sub-band are 
coded independently, and where the receiver of a link decodes 
its message on an active sub-band while treating all other 
signals on the same sub-band as interference. We assume that 
the capacity Cij(q) of link (i, j) on a sub-band q is a function 
of Xij(q), the signal-to-interference-plus-noise ratio (SINR) of 
link (i, j) over q. Denoting the transmission powers used by 
the active links on sub-band q by {P m n(q)} ( m ,n)ec q , we nave 



Xij(q) = 



2 (m,n)eC q - ■ n, , 
(m,n)^(i,j) 



G q m ,P mn (q)+N]- 



(4) 



where G* • is the path gain from m to j on sub-band q, and Nj 
is the power of the additive noise on sub-band q at j. Note that 
since the parameters {Gf^} and {Nj} are sub-band-dependent, 
this framework is particularly appropriate for networks with 
frequency selective channels and colored noise. 

Assume each node is limited by an individual power con- 
straint Pi, i.e., 



E E P *M < 



(5) 



Denote the set of power variables {Pij(q)} (i,j)ec,qeQij that 
satisfy (O by V. In Section [VI] we will design a set of power 
control algorithms that adjusts the power variables within the 
feasible region to minimize total network cost in conjunction 
with congestion control and routing. 



B. Traffic Demands, Congestion Control, and Routing 

Let the traffic demands for the network consist of a col- 
lection W of unicast sessions. Each (elastic) session w G 
W is characterized by its fixed source-destination node pair 
(0(w), D(w)) and demand rate f,„F°l We model the traffic as 
fluid flows. Assume congestion control is exercised at each 
source node. That is the source node O(vo) can control the 
rate at which w's traffic comes into the network. Denote the 
actual admitted rate of w by r^, also referred to as the end-to- 
end flow rate of w. Thus, the rate of the rejected traffic of w, 
denoted by F w , is f w — r w . After the admitted flow rates are 
determined by the congestion control at the source nodes, the 
flows entering the network are routed on (potentially) multiple 
paths from the source to destination. Let fij(w) denote the 
rate of session w traffic routed through link (i, j). The session 
flows satisfy the flow conservation constraints, i.e., 



0{w) 



(6) 



E M w ) = E /«w = : *<H. i * °W- fl ( 

jeo, keo. 



and 



f ij (w)=0, VjeOi, i = D(w). 



(7) 



(8) 



Here, U(w) denotes both the total incoming and outgoing flow 
rates of session w at an intermediate node i. If we consider 
the rejected flow F w as being routed on a virtual "overflow" 
link from 0(w) to D(w) [29], [30], the flow conservation 
constraint involving both real and virtual outgoing flows from 
the source node is given by 



jeo, 



fij(w) + F u 



(9) 



Thus, later on we will incorporate congestion control in a 
unified framework with routing. Let the set of flow vectors 
(Fw, (fi](w))(ij) e c) satisfying ©-(O be denoted by T w . 

The total flow rate on link is = J^wew /»j i w )- 
To route a flow of rate F,j from i to j, node i can split the 
traffic onto all active sub-bands of link and transmit them 
simultaneously. Let the rate of flow assigned to sub-band q £ 
Qij be Fij(q). We have = Z) 9eQy Fij(q). Hence, 



E F *M = E 

qGQij wEW 



(10) 



Note that the flow on each sub-band may consist of traffic of 
one or more sessions. However, the traffic split conserves the 
total rate of any session going through link Therefore, 

10 We assume that user w gains no extra utility by transmitting at a rate 
higher than f w . So without loss of optimality, r w can be taken as the effective 
maximal incoming rate demanded by w. 
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at the receiver end of the link, node j collects each session 
w's traffic of rate fij(w). 

One can think of the routing scheme as a two-step process. 
The first step determines inter-node routing, i.e., a feasible 
flow vector (F Wl {fij(w))uj)ec) G -Fw is found for every 
to G VV, yielding The second step amounts to 

intra-node routing: each node i routes the flow onto the 
active sub-bands q G Qij available from % to j. 



C. Network Cost and Optimal Resource Allocation 

The cost on each (real) link is made up of the costs 
Df- incurred on its active sub-band(s) q 6 Qij . Each Df- is a 
function of iy (q), the SINR of link on q, and Fij (q), the 
flow rate on the active sub-band q of link Specifically, 
define 

D!j=D(x ij (q),F i j(q)). 

Assume that D(x, F) is decreasing and convex in x for fixed 
F, and that it is increasing and convex in F for fixed x. To 
impose the hard capacity constraint F < C(x), where C(x) 
denotes the link capacity achieved by having SINR = x, 



2(x)-F 



where C(x) 



one can define D(x, F) = oo for F > C(x). An example of 
such cost functions io D(x,F) = tjt 
R\og(Kx) with constants R,K > 0@ 

We associate an increasing and concave utility function 
U w (r w ), < r w < f w , to each (elastic) session w. Thus, 
rejecting a flow of F w from the network incurs a utility loss of 
U w (f w ) — U w {f w —F w ). This utility loss can also be interpreted 
as the cost D W (F W ) = U w (f w ) — U w (r w — F w ) on a virtual 
overflow link [29] with flow rate F w . Note that by definition, 
D W (F W ) is increasing and convex in F w . 

We solve for the feasible flow and power variables that 
jointly maximize the net utility, i.e., total utility minus total 
cost on real links. With the above observation, the maximiza- 
tion of net utility is equivalent to the problem of minimizing 
the total network cost E = J2(i.j)ec D i 3 + J2 w ew D ™ 
consisting of costs on both real and overflow links. Formally, 
the Multi-Radio Minimum Cost Resource Allocation (MCRA) 
problem is 



minimize 



subject to 



J2D W (F W )+J2 E D( Xij (q),FiM) 
wew (i,j)eCqeQij 

© and (P;j(?))(,,j)e£,<?eQ« e ^ 

m and (F w , (/«(«>)) W ) 6J c) e r w , Vw £ W. 



"For modelling purposes, we assume that on the transmitter side of each 
link there is one queue for every active sub-band of that link. Specifically, the 
arrival rate to the queue of link on sub-band q is Fij(q), the aggregate 
flow routed onto through q. The cost function D(x,F) = q^_ f 

gives the expected delay in an M/M/l queue with arrival rate F and service 
capacity C(x). In our present setting, F is the sub-band specific link flow 
rate and C(x) is the sub-band specific link capacity. By the Kleinrock 
independence approximation and Jackson's Theorem, the M/M/l queue is a 
good approximation for the behavior of individual links on active sub-bands 
when the system involves Poisson stream arrivals at the entry points, a densely 
connected network, and moderate-to-heavy traffic load [29], [31]. 

12 The constants R, K correspond to the bandwidth and processing gain 
of the system. The capacity formula C(x) = Rlog(Kx) is a good 
approximation for systems with a high processing gain, e.g. CDMA networks. 



We will devise a distributed scheme that iteratively adjusts 
power control, routing and congestion control on a node-by- 
node basis so as to find the optimal solution to the MCRA 
problem. 

D. Control Variables and Optimality Conditions 

To update transmission powers, traffic routes, and traffic 
input rates, we define the following power control and routing 
variables. These optimization variables are node-based and 
local in the sense that they are independently controlled at 
individual nodes. 

1) Power Control Variables: Let Pi(q) be the total power 
used by node i on its active sub-band q. Let Pij{q) be the 
power used by node i on its outgoing link over an active 
sub-band q. Define 



constrained by 



< Pi(q) < 1, ^ L 



(11) 



Also define 



subject to constraints (let rjij(q) = for q £ Q«) 

<% (<?) <1, ^77y(<7) = l. (12) 

With the above definitions, Xij (q) can be expressed in terms 
of the power control variables as 



(13) 



njtj Vin(q) + J2m^i G^jPrnPrniq) + N] 

We use INij(q) to denote the total interference and noise 
power of link on sub-band q, i.e. the denominator of 

crj. 

Define p t = (pi(g)) ge oc ( and r] i = {r}ij{q))jeo uq eQ ir 
We now compute the partial derivative of the network cost E 
with respect to each pi{q) and rjij(q). These partial derivatives 
are useful for characterizing the optimality conditions and for 
iterative adjustment of the power control variables. We have 



dE 
driij (q) 



.n£Oi 



IN in (q) 



(14) 

where we have used shorthand notation (£)*.)' for 
dD(x in (q), F in (q)) /dx m (q) and defined 

^(1 + %(<?)) 



5 m M = {DlJ 



IN ij(l) 



(15) 



Note that Si]ij(q) as well as dE/dr/ij(q) involves only local 
measures of i. The partial derivative of E with respect to Pi(q) 
is given by 



Spi(q) 



a dE 



dpi{q) 



Pi { E g^msgI + J2 SmMrniq) 



(16) 



s 



where 



MSGl 4 J2 (Dl 
meO n 



iN mn ( q y 



(17) 



is the message that node n needs to send to any other node, 
say i, which is active on q in order for i to compute 6pi(q) 
via ( TToT l. Specifically, the message exchange works as follows. 

Power Control Message Exchange Protocol: Let each node 
n keep an up-to-date measure MSG^ for each sub-band q, 
where MSG^ is derived by assembling the measures 

/ p.q y ~GmnPmn(q) _ r p. q y ~ ^mwCg) 



IN mn (q) 2 m " GmnPmnil) 

on all its active incoming links (m, n) on g, and summing them 
up (cf. (fTT])). Note that MSG^ is nonnegative and it is zero 
if n has no active incoming link on q. If MSG^ > 0, node n 
broadcasts it to the whole network via a flooding protocol. This 
control message generating process is illustrated by Figure [3] 
If q E OCi, node i collects MSG^ and processes it according 



Dt 




^ - ( * 
MSG" = V D q ' X " mKq) 

>„„, G « p , , 



Fig. 3. Power control message generation. 



to the following rule. Node i multiplies MSG^ with the path 
gain G| n O It further adds the product to the value of local 
measure 5r)i n (q) ■ T]i n (q) if the origin n of the message is a 
neighbor of i and (i, n) is active on sub-band q. Finally, node 
i adds up all the processed messages, and this sum multiplied 
by Pi equals Spi(q). Note that this protocol requires only one 
message from each node on each sub-band Moreover in 
practice, node i can ignore the messages generated by distant 
nodes, because they contribute very little to 5pi(q) due to the 
negligible multiplicative factor G\ n on MSG^ when i and n 
are far apart (cf. (fToll). 

2) Routing Variables: Routing variables were first intro- 
duced by Gallager [32] for wireline network routing problems. 
Here, we define routing variables in a similar fashion. In 
addition to inter-node routing, however, routing variables here 
also perform the function of congestion control and intra-node 
routing. 

Recall that congestion control is equivalent to routing a 
portion of traffic demand on a virtual overflow link directly 

13 In a symmetric duplex channel, ~ G'L and node i may use its own 
measurement of G ni in the place of G iri . Otherwise, it will need channel 
feedback from node n to calculate G| . 

14 To be more precise, only the nodes having at least one active incoming 
link on a certain sub-band need to provide a message for that sub-band. 



from the source to the destination. Let i be the source node 
of session w, define overflow routing variable 

I A A W 
<Pw = — , 

which is constrained by <<f> w < 1. The overflow rate is then 
controlled by <fi w as F w — f w <fr w , and the end-to-end flow rate 
is given by r w — f w (l — <j> w ). Routing variables associated 
with a real link are defined by 

which gives the fraction of the incoming flow of session w at 
node i that is routed onto link At any node i except the 

destination, fcj (w) of all j E Oi satisfy 

< (f>ij{w) < 1 and ^ <t>ij{w) = 1. (18) 

It is easy to see that the routing variables tp w , (w) of all 
E C and w E W uniquely determine the inter-node flow 
patterns (F w , {fij{w))(i,j)ec) of all sessions w, and hence the 
total flow rate on links Fij = 2~2 W fij( w )- Now let the flow 
allocation on active sub-bands be specified by the intra-node 
routing variables defined by 

■fij 

For any link we must have 

0</i y (g)<l and ^ lHj (q) = l. (19) 
qeQ i3 

The link flow rate on an active sub-band is therefore given by 

Fijiq) = lHj{q) ^2 ti(w)(j)ij(w), q E Qij. 

For notational brevity, denote by 4>i(w) the vector 
(i>ij(w))jeOi if i ^ 0{w),D(w), or {(j)w,{ ( t>ij{' w ))jeO i ) if 
i = 0(w). Also denote by /^^ the vector (p*ij(q))j£O i ,qeQ ij - 
We illustrate the use of the above routing variables by looking 
at a source node i and one of its outgoing links in Figures 
H and [5] In Figure 2J the overflow link is marked by a solid 
line with hollow arrow. In Figure [5] we assume that link 




Fig. 4. Overflow and inter-node routing at a source node. 

is active on two sub-bands qi and q2, and that the inter-node 
flow rate from i to j is Fij . 

We now compute the partial derivative of E with respect to 
those routing variables. First we have the derivative of E with 
respect to the overflow routing variables 



dE 



DL - 



dE 

dt„, 
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Fig. 5. Intra-node routing on link 



where J-p is the marginal cost of increasing the end-to-end 
flow rate of session w with all other variables held constant. 



dE 
dti(w) 



which is the 



This is a special case (i = 0(w)) of 
marginal cost of increasing the incoming flow rate of w at node 
i while keeping all other variables constant. These marginal 
costs are computed recursively as in [32]: 

dE 



&ti(w) 
and for i ^ D(w), 



0, if % = D{w), 



dti(w 



E 0ii(«>)$(fe(lD 



dtj(w) 



(20) 



where we have used shorthand notation D% ' for 

dD(xij(q), Fij(q))/dFij(q) and defined marginal routing 
cost indicators 

dE 



;H = E /'•.•' : '/ ,/; '/--,/ 



dtj(w) 



(21) 



Using 6<t>ij(w), we can easily write out the partial derivative 

of E in (f>ij (w) as 

OE 

7T7 — 7 — r = t l (w)d(j) lj (w). (22) 
d<t>ij(w) 

Finally, the partial derivative of U in /ijj (g) is given by 

dE 



FD q ' 



(23) 



Notice that q^ E ^ depends on local measures only. However, 
dE and a , dE , are both tied to the value of -^fr of i's 

0{p^(w) Otj{W) 



downstrearro neighbors j, which in turn depends on 
of j's downstream neighbors k. Thus, we need a sequential 
message passing of the marginal routing costs, from the 
destination upstream to the source, to permit every node to 
acquire the partial derivatives in its local routing variables. 

Routing Message Exchange Protocol: In [32], the rules 
for propagating the marginal routing cost information are 
specified. In order for node i to evaluate the terms 8<f>ij (w) in 
(|2T1 >. it needs to collect local measures D q F .. from all q G Qij 

15 Given a routing configuration {<j>ij( w )}(i,j)£C °f session w, node j is 
said to be downstream to i if there exists a path (ji,j2), • •', (jn,j) 

such that <pij 1 , <feii' 2 > ' ' ' ■ fynj are a " positive. We say i is upstream to j 
with respect to session w if j is downstream to i. 



dE 
dt k (u 



as well as reports of marginal costs dD/drj(w) from its next- 
hop neighbors j G Oi, for all traversing sessions w G W. 
Moreover, it is responsible for calculating its own measure of 
marginal cost d f® w ^ with respect to every session w according 
to d20l i. and then providing the measure to its upstream neigh- 
bors with respect to the session w. The sequential message 
passing terminates if and only if the routing pattern of the 
session contains no loops, which is guaranteed by the blocked- 
node-set technique developed in [32], [33]. 

3) Conditions for Optimality: The power and routing con- 
figuration that solves the MCRA problem can be characterized 
in terms of the marginal power and routing costs as follows. 

Theorem 1: For an instance of optimization variables 
{p i ,iliA'Pi(' w )}w&W,IJ'i}ieAr to be optimal, the following 
conditions are necessary: For all i G N and w G W with 

U(w) > 0, 



dE 



mm 

jeOt 



dti(w) 
Moreover, for all w G W, 



mm 



Ft 



dE 
dtj(w) 




= D' W , 
> D' 



if 4> w = 0, 
if < <j) w < 1, 
if 6,„ = 1. 



For alH G J\f and q G OCi, there exists a such that 



= A, 

> A,; 



if Pi (q) > 0, 
if pi(q) = 0, 



(24) 



(25) 



(26) 



and the constant Aj satisfies 

{ > 0, if E g ee>c» Piil) 
Xil =0, if 0<J2 qeOCzPt (q) < 1, (27) 

I < 0, if EgeOCi Pi(«) = L 

Furthermore, for q G OCi such that pi{q) > 0, there exists a 
constant ~fi(q) such that for all j G with q G <2y, 



$Vij{q) 



(28) 



= 7i(?)) if liife) > 0. 

>7i(g), ifj7y(g) = 0. 

We sketch the proof of Theorem Q] in Appendix [C] First 
order conditions as in (l24li-(|2"8Tl are in general only necessary 
for a configuration to be optimal. However, they are further 
sufficient if the link cost function D(x, F) has certain con- 
vexity properties as we discuss next. 

Let x generically represent the SINR of a link on one of 
its active sub-bands. By x is a function of the power 
variables of all active links on the same sub-band, e.g. 
Xij(q) = h^((P mn (q))( m ^) E c q )' Thus, D\- is a function of 
the vector P(q) = {P m n{q))( m ,n)ec q and F tJ (q). It turns 
out that a characterization of the sufficient conditions for 
optimality requires D 9 j to be jointly convex in P(q) and 
Fij(q). However, it can be shown that such a property is 
impossible given the assumption that Df- is decreasing in 
Xij(q). 

One way to remedy this problem is to work with log-power 
variables, first introduced in [1], [2]. For each (m, n) G C and 
q £ Qmn, define S mn (q) = lnP mn (q). Consider as a 
function of S(q) = {S mn {q))( m ,n)ec q and F i 3 {q)- If D % is 
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jointly convex in S(q) and F^j (q), the conditions in TheoremQ] 
are enough to imply the optimality of a configuration. The joint 
convexity property holds if and only if the link cost function 
D(x, F) satisfies the following. 

Lemma 6: For all G C and q G Qij, Df- is jointly 

convex in S(q) and Fij(q) if and only if the matrix 



M 



D'^x 2 + 
D" 



D'x D'l 



n" 



(29) 



is positive semidefinite for all (x,F). 



In (|29l , D' x and Z?" denote the first and second partial 
derivatives of D(x, F) with respect to x, D' F denotes the 
second partial derivative of D with respect to F, and D'£ F 
denotes d 2 D / dxdF . It can be shown that the cost function 
D(x,F) — C (*}_ F with C(x) = R\og(Kx) satisfies the 
condition in Lemma [6] Limited by the available space, we 
skip the proof of Lemma [6] We prove in Appendix [D] that if 
the condition in Lemma [6] holds, the conditions stated in the 
next theorem are sufficient for optimality. 



Theorem 2: If the matrix M in d29l of link cost function 
D(x,F) is positive semidefinite, then conditions (I24li-(|28T> are 
sufficient for {p i , rj i , {</>i(w)} we y\>, PijieN t0 be optimal if 
(l24l holds at every node i ^ D(w) whether ti(w) > or 
not, and if (f28j holds for all i G Af and q G OC t whether 
Pi(q) > or not. 



VI. Node-Based Multi-Radio Power Control and 
Routing Algorithms 

We develop a set of scaled gradient projection algo- 
rithms [34] by which individual nodes adjust their local 
power control and routing variables iteratively to achieve a 
global configuration satisfying the optimality conditions in 
Theorem |2] Each algorithm at a node updates the appropriate 
vector of local optimization variables, e.g. p { , 77^ 4>i{ w ) or 
p t , such that the updated vector results in a lower cost with 
all other variables held constant. This is achieved by updating 
the old vector in the opposite gradient direction scaled by 
an appropriate positive definite matrix, and projecting the new 
vector back into the feasible set whenever it falls outside. This 
technique has been applied to, for example, optimal routing in 
wireline networks [33] and optimal power control and routing 
in single -radio wireless networks [3]. 

A. Power Control Algorithms 

We develop two sets of power control algorithms which 
let each node i iteratively adjust the vector (T)ij(<l))j&Oi an d 
the vector (pi{q))q£OCi- First notice that if q G OCi but 
L q contains only one of i's outgoing links, say i]ij(q) 
must be equal to 1, and hence is not variable. Otherwise, if 
Vi(l) — {Vij (l))(i,j)&c q has more than one positive element, 
it is updated by the following scaled gradient projection 
algorithm: 



2? (?) 



(30) 



Here, the superscripts k, k + 1 are the iteration indices, Qi(q) 
is a positive definite scaling matrix, Srj^q) is the vector 
(8rjij(q))(ij) e c q , and Hqfcfa) is the projection operation onto 
the feasible set of rj i (q) (cf. (U2l ) relative to the norm induced 
by Q* (q)0To implement the algorithm (l30l l. node i needs the 
current value of Srj i (q), which by (fl5T l is easily computable 
from local measures. 

The algorithm used by node i to iteratively update p i is as 
follows: 



(31) 



It has almost the same form as d30T > except that the feasible 
set for the projection is defined by (fTTl i. However, notice 
that each component of bp\ = (8pi(q) k ) q <£OCi depends on 
measures from all active links on sub-band q. Therefore, prior 
to implementing (|3T1 >. node i needs to collect the appropriate 
power control messages to determine Spf , as described by the 
message exchange protocol in Section IV-D.ll 

B. Routing Algorithms 

We now present the algorithm used by a node i to update 
its inter-node routing vector 4>ii w ) f° r a traversing session w. 
An iteration of the inter-node routing algorithm has the form 



d> k+ \w) := k fc H - (Atf H)- 1 • 54> k {w) 



+ 



■ (32) 

Here, 5<fi k (w) = (50«('w))jeOj is the vector of current 
marginal routing cost indicators and M k (w) is a positive 
definite matrix used to scale the descent direction. The fea- 
sible set associated with the projection operation in (l32l l is 
prescribed by the constraint ([LSl i and an additional requirement 
that (f>ij(w) = for all j G B k (w), where B k {w) is the 
blocked node set of node i relative to session w. The meaning 
and use of the blocked node set were briefly discussed in 
the routing message exchange protocol in Section IV-D.2I The 
exchange protocol also gives every node enough information 
to calculate the current 84>i{w) prior to each iteration of ( l32l . 

For the source node 0(w) of a session w, the algorithm 
(l32t is applied to the vector of routing variables associated 
with real outgoing links. The overflow routing variable 4>w is 
updated by 

+ 



dL 



8E 
dr w 



(33) 



< 1. 



for which the projection is onto the feasible set < 
The gradient is given by subtracting J^-, which is computable 
after the routing message exchange, from the local measure 

K- 

Finally, we come to the algorithm that node i uses 
to iteratively adjust the intra-node routing vector p,^ = 
(Viji^qeOCi applied to an outgoing link The intra- 

node routing update is iterated as 



fc+i 



(T k y 



(34) 



16 



In general, [x]~f, = arg min^^yr^x — x) 1 ■ M ■ (x — x), where T is 



the feasible set of x. 
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Here, dDp i:j is the vector of partial derivatives dD^j/dFij(q) 
of all q G Qij, which are purely local measures. The scaling 
matrix T* is positive definite. The feasible set that the 
projection refers to is defined by ( fT9b . 



C. Convergence of Algorithms 

Having described the node-based power control and routing 
algorithms, we now state the main convergence result in the 
following theorem. 

Theorem 3: Given a feasible spectrum allocation {C q } qG Q, 
let {»7°}, {Pi}, {(fi^w)}, {Hi} be any feasible initial trans- 
mission power and routing configuration with finite cost. 
Then with appropriate scaling matrices, the update sequences 

R fc }r=i> {P"}k=v {<t>lH}?=i and {tf}k=i generated 
by the algorithms (l30i >. PTT l. (l32l and d34t converge, i.e., 

Vi -» nt. Pi -» Pi <t> k i{w) -» 4>*H, and p% -> p* 
for all i G M and w G W as k — > oo. Furthermore, the 
limiting configuration {r/*}, {p*}, {<p*(w)}, {p*} satisfies 
the conditions in Theorem [2] and is a jointly optimal solution 
to the MCRA problem if the link cost function D(x, F) 
satisfies the condition in Lemma [6] 

The proof hinges on the fact that, by using appropriate scal- 
ing matrices, every iteration of the algorithms ( f30b . (l3Tl l. (l3~TT l 
and ( f32l reduces the total cost of the MCRA problem until the 
optimality conditions in Theorem [2] are achieved. Finding the 
appropriate scaling matrices, however, is a major challenge. 
One approach to this problem is to choose the scaling matrices 
so that they upper bound the Hessian matrices with respect 
to the updated variables. In this way, the algorithms closely 
approximate Newton's method, hence enjoying fast rate of 
convergence while simultaneously guaranteeing convergence 
from all initial conditions. This method has been successfully 
adopted in the power control and routing algorithms for single- 
radio wireless networks in [3], and can be generalized to the 
present context. Due to the limited space, however, we skip 
the details. 

It is worth noting that convergence does not depend on any 
particular order of running the algorithms at different nodes. At 
any time, any node can update any set of its local optimization 
variables via the corresponding algorithm. All that is required 
for convergence is that each node iterates every algorithm 
until the adjusted variables have marginal costs satisfying 
conditions in Theorem [2] Finally, we note that convergence 
occurs from any initial configuration with finite cost. These 
features are crucial to the applicability of these algorithms 
to large networks which lack the ability of scheduling and 
synchronizing node operations. 

Extensive simulations indicate that our algorithms are adap- 
tive to time-varying network state, including channel fading, 
network topology, and traffic demand. Because every iteration 
of any algorithm always reduces the total cost under the 
current network condition, our scheme is able to constantly 
readjust routing and transmission powers towards the optimum 
that slowly shifts over time due to the network change. 



VII. Conclusion 

We have developed an integrated cross-layer resource allo- 
cation scheme for general wireless multi-hop networks. To 
satisfy the fundamental duplexing constraints, our scheme 
first finds a feasible spectrum allocation by (1) dividing the 
whole spectrum into multiple sub-bands and (2) activating 
conflict-free links on each sub-band. Compared with tradi- 
tional scheduling in time, the spectrum allocation technique 
has a number of advantages in operational simplicity and 
amenability to distributed and asynchronous implementation. 
By studying an equivalent combinatorial link-coloring prob- 
lem, we found that the minimum number of sub-bands re- 
quired by a feasible spectrum allocation is given by a simple 
function of the chromatic number of the network connectiv- 
ity graph. The minimum number grows asymptotically at a 
logarithmic rate with the chromatic number, attesting to the 
good scalability of the spectrum allocation technique and its 
robustness to network topology changes. We designed a simple 
distributed and asynchronous algorithm by which a feasible 
spectrum allocation can be constructed given enough sub- 
bands. 

Given a feasible spectrum allocation, we developed an 
analytical framework and a set of node-based distributed 
algorithms for optimally allocating transmission powers and 
traffic rates on active links. Such a framework is especially 
suitable for the design of wireless networks with frequency 
selective channels. We provided the conditions that an optimal 
power control and routing configuration must satisfy. We 
then designed a set of distributed power control and routing 
algorithms using the scaled gradient projection method. These 
algorithms can be iterated at individual nodes with little 
control overhead. Finally, we demonstrated that the algorithms 
asymptotically achieve the optimal configuration regardless 
of the initial condition and the order of iterating different 
algorithms. 

Appendix 

A. Proof of Lemma |2] 

We provide the proof for the case g > [Q/2\. The other 
case g < [Q/2\ can be seen as a corollary by taking 
complements of all subsets involved in the first case. 

First notice that since g > [Q/2\, there are K = ( _ x ) 
distinct subsets C%, ■ ■ ■ , Ck of cardinality g — 1, where K > 
k = (^), the number of distinct subsets of cardinality g. 
Hence, the claim makes intuitive sense. 

Consider the bipartite graph consisting of {Bi} and {Cj} 
where a pair of Bi and Cj are connected if and only if Cj C 
Bi. In this case, Cj is said to be a child of Bi, and Bi is said to 
be a parent of Cj . It is easy to see that every Bi has g children 
and every Cj has Q — g + 1 parents. Now the claim in the 
lemma is equivalent to the existence of a complete matching 
(of size fc) of {Bi} and {Cj} in the bipartite graph specified 
above. Because g > [Q/2\ implies that g > the degree 
of any Bi, which is g, is greater than or equal to Q — g + 1, 
the degree of any Cj. In this case, a complete matching must 
exist by a corollary of Hall's theorem (cf. Corollary 13.4 of 
[35]). 
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B. Proof of Lemma \3\ 

Suppose {Bi]i e ^f is an optimal solution and HJi-^il = Q- 
If {Bi} does not have the property in the lemma, we can 
always modify it to {B*} which satisfies the property with 
Q* = ||L B* | being less than or equal to Q. Hence, the lemma 
follows. Suppose we have a feasible solution B consisting of 
a collection of N subsets {Bi)f =l . We show that based on 
{Bi}, we can construct B = {Bi}f =1 such that \Bi\ = [Q/2\ 
for all % = 1, ■ • • ,N, where Q = | U B;|. 

Define m = min; \Bi \ and M = maxj |Bj|. If m < M, then 
to < LQ/2J or M > LQ/2J or both. If to < [Q/2\, define 
M m = {i € TV : \Bi\ = m}, replace each Bj, j € TV m , by 
subset B^ with cardinality to+1 such that Bj C Bj for all j G 
M m , and Bj ^ B^, if and only if Bj ^ Bk- Such a replacement 
is possible by Lemma [2] If M > [Q/2\, define Mm = £ 
TV : |Bj| = M}, replace each Bj, j E Mai, by subset Bj 
with cardinality A/ - 1 such that Bj C Bj for all j E TVm, 
and Bj ^ B' k if and only if Bj ^ B^. Such a replacement 
can be found also by Lemma [2] Denote by TV' the subset 
of nodes whose Bi is changed (either expanded or reduced). 
Because to < M, M' is always non-empty. It can be verified 
that the new collection of subsets {{Bj}j e ^/> , {Bij^^'l is 
another optimal solution of (fJJ. If the minimum and maximum 
cardinalities of the new collection of subsets are equal, we are 
done. Otherwise repeat the above procedure until we obtain 
an even collection of subsets. The iterations terminate in a 
finite number of steps since each iteration strictly reduces the 
difference between M and m. 

C. Proof of Theorem [7] 

We show that whenever one of conditions (|24|i-(|2"8Ti is 
violated, the present configuration can be improved upon. We 
take condition (124-b as an example. Arguments for the other 
conditions are similar. By (ETT i and (EOt . we have 

dE 



and 



l e ^ J dtj(w 



dE ^ . 
— — ; — r > mm < mm 

oUlw) jeOi [qGQ 



Thus, condition d24"l i holds if and only if for all j S O 



such that 4>ij{w) > 0, 54>ij(w) < min ge g i( . 



D'l 



dE 
dt k {w) 



for all k £ Oi. Suppose condition (1241) is violated at i with 
ti(w) > 0, i.e., there exists j E Oi such that 4>ij{w) > 0, 



i(w) > min 9eaifc 



dE 
dt k (w 



a for some k E Oi 



If Fik > and 5<pik{w) > a, there must exist q E Qik with 
A*i/c(<z) > but D F ' > Dp, k ' for some other v E Qik- 
Hence, the cost of the configuration can be further reduced 
by shifting a tiny portion from ^ik(q) to fiik(v) as suggested 
by (|23T ). If 8<pik{w) = a, then 54>ij(w) > 8<pik{w), we can 
reduce the cost by shifting a tiny portion from 4*ij{w) to 
4>ik{w) in light of (1221 l, If Fik = 0, then we can always make 
54>ik(w) — a by setting all fiik{q) equal to zero but /J,ik{v) 
equal to one for one v achieving the minimum in a. Notice 
that reconfiguring {^ik{q)}qeQ ik does not change the value of 
er since by assumption Fik = 0. Then, as before, the total cost 



can be further reduced by shifting a tiny portion from <f>ij (w) 

to (f>ik(w). 

D. Proof of Theorem [2] 

Suppose {pi, t]i, {0 i (w)} u)e w, Mi} satisfies the conditions 
in Theorem [2] Suppose the configuration yields log-power 
variables {S(q)} qe Q and link flow variables (F w ) w< zw, 
( F ij(l))(i,j)eC, qeQij- Let {p u f}^ {4>i(w)}weW, £J be an- 
other feasible configuration which yields {S^q)}, (F w ) we w, 
{Fi 3 {q))(i ,j)es, qeQij - Recall that each By is jointly convex 
in S(q) = (S mn (q))( m ,n)ec q and Fijio) whi l e each D w 
is convex in F w . Moreover, the feasible sets of log-power 
variables and flow variables are both convex. Therefore, the 
cost difference under the two configurations can be bounded 

by 

£a»+ E E s&-£a.- E E D % 
> D '™( p ™ - F -) + E E D h/(^(i) - km) 

dD q - 

13 -{S mn {q)- S mn {q))- (35) 



mew 

+ E E E 

(i,j)SC q&Qij (m,n)eC q 



dSmn(q) 



We will show that the RHS is nonnegative. We can re-write 
and lower bound the first two summations by the series of 
equalities and inequalities on the top of the next page. 

For equality (a), we used the relation that 

weW wy ' (i,j)e£ge&y 

and appended terms [tj(w) - T, i7 tD(w) U(w)^ij(w) 

which are all equal to zero by flow conservation constraints. 
Equality (b) is obtained by rewriting r w and f w as f w (l — 
4> w ) and r w (l — <p w ), respectively, and reorganizing the terms. 
Condition d24"l i implies that 

dE \ ^ \ - _ dE 

dun - \ \ n > ■ JJU^r 

by which we obtain (c). The final inequality follows from 
condition (|25T >. 

Next we lower bound the third summation on the RHS of 
( |35"1 >. To begin with, we switch the roles of (m,n) and 
in the summation to rewrite it as 
fjjjq 

E E i^iSiM-SiM)- 



The partial derivative is computed as 



dD q mn 

dS l3 (q) 



if = (m,n), 



- D XmJ X mn(q) G IN^{q) ' ° therwise - 



Thus, we can expand the summation as shown in the second 
block on the next page. 

For the first summation on the RHS of d36l l, we use the 
convention that 0/0 = 1 and y/0 > 1 for all y > 0. As 
a consequence, the summand vanishes if fjij(q) — = 
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(a) 



(<>) 



> 



D' Fm ■ (F w - F w ) + E D k 3 ' ■ (?) - (g)) 



E ^ ■ »■ 



E E E [^ i /-ti(»)^MAy(?)] 

mew [(i,i)e£?eey 



dE 
dt w (w) 



' t w (w) 



E E 



dE 
dtj(w) 



D' w <j> w + (1 — 0™) 



9£ 



+ E E **(«o 

u)£W I i^O(w), D(w) 



E <^H E D F i} 'P*i(<l) 
\q€Qij 



dE 



dE 



E> 
E> 



D w (j>w + (1 — 0id) 

Ok 



jeOi 
3£ 



dt w (w) 
dE 



d' w 4>w + (i — 



dtj(w) j dU(w) 
dE 



dt w (w) 



+ E E &(«) 

tu£W I i^O(w),D(w) 



dE dE 



dti(w) dU(w)_ 



dt w (w) 



> 0. 



E E E 



''^(SiM-SiM) 



dS tj (q) 



E E 

(i,j)€C qeQij 



E E 
E E 

(»,j')e£«eey 



e: -jl'^w % p<i ff + ^- ■ v ) 



(m,n)e£ g 
.(m,n)/(i,j) 



IN mn (q) 



In 



E -^/^(^ i^i + £>g i /:r ij ( 7 n I +.,-,,.(,/)) 



(m,n)e£g 



/iV m „(g) 



In 



Py(«)In 



£A E *&) E 



E G in. MSG n + 



gj(g) 



(36) 



0. Moreover, for those € £ g but r)ij(q) = 0, we can Now consider the second summation on the RHS of d36b . 

lower bound the summand by replacing 8r)ij(q) with ji(q) For each i € J\f, let the summands with pi(q) = be lower 
(cf. condition d28l >) to get bounded by 

' r„\\ i ft (?) 



^ By condition ( l26*l >, summands with /9j (g) > are equal to 



It is implicit from ( [Tol l and d28l ) that Spi(q) = pi{q) 
Pi [Ene./V G\ n MSGl + ji (g)]. We thus can lower bound the ft (?) A * ln ■ 

first summation in (|36*1 i by 

Therefore, the whole summation is lower bounded by 

> E E ft(^ft(«) E ^^)f!44-i)=o. > E E ft(^Y^-i 



(*>j')eAr / 
The inequality follows from the relation ln x < x — 1 and that = \^ \ { I \^ Pi(g) — Pi (g) > 0. 

Pi(q)8pi(q) < due to conditions ( f26b and (f27b . ieA^ VgeoCi rjeoCi 
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The second to last inequality is obtained by noting that 
Pi(q)\ < due to condition d27| ) and the identity In a; < x—1. 
The last inequality follows also from condition d27| ). 

Thus, we have shown that the cost incurred by an arbitrary 
feasible configuration is greater than or equal to the cost by 
a configuration that satisfies the conditions in Theorem [2] So 
the proof is complete. 
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